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Abstract 

The simple assembly line balancing problem (SALBP) concerns the assignment of tasks 
with pre-defined processing times to work stations that arc arranged in a line. Hereby, 
precedence constraints between the tasks must be respected. The optimization goal of the 
SALBP-2 version of the problem concerns the minimization of the so-called cycle time, 
that is, the time in which the tasks of each work station must be completed. 

In this work we propose to tackle this problem with an iterative search method based 
on beam search. The proposed algorithm is able to obtain optimal, respectively best- 
known, solutions in 283 out of 302 test cases. Moreover, for 9 further test cases the 
algorithm is able to produce new best-known solutions. These numbers indicate that the 
proposed iterative beam search algorithm is currently a state-of-the-art method for the 
SALBP-2. 

1 Introduction 

The class of problems known as assembly line balancing problems (ALBPs) concerns the 
optimization of processes related to the manufacturing of products via assembly lines. Their 
importance in the industrial world is shown by the fact that much research efforts have been 
dedicated to many different types of ALBPs during the past 50-60 years ^S^iSj. The specific 
problem considered in this paper is the so-called simple assembly line balancing problem 
(SALBP) [21], a well-studied scientific test case. An assembly line is composed of a set of 
work stations arranged in a line, and by a transport system which moves the product to be 
manufactured along the line. The product is manufactured by executing a given set of tasks. 
Each of these tasks has a pre-defined processing time. In order to obtain a solution to a given 
SALBP instance, all tasks must be assigned to work stations subject to precedence constraints 
between the tasks. In the context of the SALBP, all work stations are considered to be of 
equal size. Moreover, the assembly line is assumed to move in constant speed. This implies a 
maximum of C time units — the so-called cycle time — for processing the tasks assigned to each 
work station. The SALBP has been tackled with several objective functions among which the 
following ones are the most studied ones in the literature: 
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• Given a fixed cycle time C, the optimization goal consists in minimizing the number of 
necessary work stations. This version of the problem is refered to as SALBP-1. 

• Given a fixed number m of work stations, the goal is to minimize the cycle time C. The 
literature knows this second problem version as SALBP-2. 

The feasibility problem SALBP-F arises when both a cycle time C and a number of work 
stations m is given and the goal is to find a feasible solution respecting C and m. In this 
work we will deal with the SALBP-2 version of the problem. 

For what concerns the comparison between SALBP-1 and SALBP-2, much of the sci- 
entific work has been dedicated to the SALBP-1. However, also for the SALBP-2 exists a 
considerable body of research papers. An excellent survey was provided by [21]. Approaches 
for the SALBP-2 can basically be classified as either iterative solution approaches or direct 
solution approaches. Iterative approaches tackle the problem by iteratively solving a series of 
SALBP-F problems that are obtained by fixing the cycle time. This process is started with a 
cycle time that is set to some calculated upper bound. This cycle time is then decremented 
during the iterative process, which stops as soon as no solution for the corresponding SALBP- 
F problem can be found. In contrast to these indirect approaches, direct approaches intend 
to solve a given SALBP-2 instance directly. 

Heuristic as well as complete approaches have been devised for the SALBP-2. Among the 
existing complete methods we find iterative approaches such as the ones proposed in [9l ED] 
but also direct approaches such as the ones described in |19]. I12j. Moreover, the performance 
of different integer programming formulation of the SALBP-2 has been evaluated in [16j . The 
currently best-performing exact method is SALOME-2 [12]. Surprisingly this exact method 
even outperforms the existing heuristic and metaheuristic approaches for the SALBP-2. The 
most successful metaheuristic approach to date is a tabu search method proposed in [22]. 
Another tabu search proposal can be found in [6J. Other metaheuristic approaches include 
evolutionary algorithms [21 [25l [H] and simulated annealing [10] . Finally, a two-phase heuristic 
based on linear programming can be found in [23], whereas a heuristic based on Petri nets 
was proposed in [TT] . 

Contribution of this work. Subsequently we propose to tackle the SALBP-2 by means 
of an iterative approach based on beam search, which is an incomplete variant of branch 
& bound The resulting iterative beam search algorithm is inspired by one of the current 
state-of-the-art methods for the SALBP-1, namely Beam- AGO [3]. Beam- AGO is a hybrid 
approach that is obtained by combining the metaheuristic ant colony optimization with beam 
search. In this work we propose to use the beam search component of Beam-AGO in an 
iterative way for obtaining good SALBP-2 solutions. Our computational results show indeed 
that the proposed algorithm is currently a state-of-the-art method for the SALBP-2. It is able 
to obtain optimal, respectively best-known, solutions in 283 out of 302 test cases. Moreover, 
in further 9 cases the algorithm is able to produce new best-known solutions. 

Organization of the paper. In Section [2] we present a technical description of the tackled 
problem. Furthermore, in Section [3] the proposed algorithm is described. Finally, in Section |4] 

^The interested reader may note that we proposed a very similar approach for a more general problem, the 
assembly line worker assignment and balancing problem (ALWABP) [S]. 
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we present a detailed experimental evaluation and in Section [5] we conclude our work and 
offer an outlook to future work. 



2 The SALBP-2 

The SALBP-2 can technically be described as follows. An instance (T, G, m) consists of three 
components. T = {1, . . . , n} is a set of n tasks. Each task i €T has a pre-defined processing 
time ti > 0. Moreover, given is a precedence graph G = {T,A), which is a directed, acyclic 
graph with T as node set. Finally, m is the pre-defined number of work stations which are 
ordered from 1 to m. An arc lij € A indicates that i (z T must be processed before j € T. 
Given a task j € T, Pj C T denotes the set of tasks that must be processed before j. A 
feasible solution is obtained by assigning each task to exactly one work station such that 
the precedence constraints between the tasks are satisfied. The objective function consists in 
minimizing the so-called cycle time. The SALBP-2 can be expressed in the following way as 
an integer programming (IP) problem. 



subject to: 



min z (1) 

m 

J2^^s = 1 Vier (2) 

s=l 

m m 

sxis < ^ sxjs Vj eT,i e Pj (3) 

s=l s=l 

'^tiXis < z s = l,...,m (4) 

Xis € {0,1} Vi € r, s = 1, . . . , m (5) 

z > (6) 

This IP model makes use of the following variables and constants: binary variable 

which is set to 1 if and only if task z G T is assigned to work station 1 < s < m. The 
objective function (1) minimizes the cycle time z > 0| The constraints (2) ensure that 
each task i € T is assigned to a single work station 1 < s < m. Constraints (3) reflect the 
precedence relationships between the tasks. More specifically, if task j € T is assigned to a 
work station 1 < s < m, all tasks i € Pj must be assigned to work stations 1 < s' < m with 
s' < s. The constraints (4) ensure that the sum of the processing times of the tasks assigned 
to a work station 1 < s < m do not exceed the cycle time z. Note that this model is not 
necessarily the most efficient IP model for solving the SALBP-2. An evaluation of several 
different models can be found in [16j . 



2.1 Solution Representation 

The following solution representation is used for the description of the algorithm as given in 
Section[3l A solution S is an ordered list S = {Si, . . . , Sm) of m sets of tasks, where Si denotes 

■^Note that we refer to the variable cycle time of the IP model as z, while fixed cycle times are denoted by 

C. 
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the set of tasks that are assigned to the i-th work station. Abusing notation we henceforth 
call Si a work station. Note that for a solution S to be valid the following conditions must 
be fulfilled: 

1. UI^i Si = T = {1, . . . ,n} and CiiLi = 0, that is, each task is assigned to exactly one 
work station. 

2. For each task j G Si it is must hold that Pj C \J\^i S^- This ensures that the precedence 
constraints between the tasks are not violated. 

2.2 Reverse Problem Instances 

The reverse problem instance (T, Cjm) with respect to an original instance {T,G,m) is 
obtained by inverting the direction of all arcs of G. It is well-known from the literature [5T] 
that tackling the reverse problem instance may lead an exact algorithm faster to an optimal 
solution, respectively, may provide a better heuristic solution when tackled with the same 
heuristic as the original problem instance. Moreover, a solution 5'" to the reverse problem 
instance (T, C" , m) can easily be converted into a solution S to the original problem instance 
(T, G, m) as follows: 

Si := Sl^_i^i for z = 1, . . . , m (7) 

3 Iterative Beam Search 

As mentioned in the introduction, the basic component of our algorithm for the SALBP-2 
consists of beam search (BS), which is an incomplete derivative of branch &: bound [15j. 
Initially BS has especially been used in the context of scheduling problems (see, for exam- 
ple, [TTl El [24] ) . To date only very few applications to other types of problems exist (see, for 
example, |13[ [Hll]). In the following we briefly describe how one of the standard versions of 
BS works. The crucial aspect of BS is the parallel extension of partial solutions in several 
ways. At all times, the algorithm keeps a set B of at most /c^w partial solutions, where B is 
the so-called beam, and /Cbw is known as the beam width. At each step, at most fc^xt feasible 
extensions of each partial solution in B are selected on the basis of greedy information. In 
general, this selection is done deterministically. At the end of each step, the algorithm creates 
a new beam B by choosing up to /Cbw partial solutions from the set of selected feasible ex- 
tensions. For that purpose, BS algorithms determine — in the case of minimization — a lower 
bound value for each extension. Only the maximally fcbw best extensions — with respect to 
these lower bound values — are included in B. Finally, if any complete solution was generated, 
the algorithm returns the best of those. Note that the underlying constructive heuristic that 
defines feasible extensions of partial solutions and the lower bound function for evaluating 
partial solutions are crucial for the working of BS. 

In the following we first present a description of the implementation of the BS component, 
before we describe the algorithmic scheme in which this BS component is used. 

3.1 The Beam Search Component 

The BS component described in this section — see Algorithm [1] for the pseudo-code — is the 
main component of the proposed algorithm for the SALBP-2. The algorithm requires a 
problem instance (T,G,m), a fixed cycle time C, a beam width /c^w, and a maximal number 
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Algorithm 1 Beam search 



input: an instance {T,G,m), a fixed cycle time C, a beam width /c^w, and fc^^^t 
/ := 

Initialization of an empty solution S 
B := {S] 

-Bcompl •= 

while 5 7^ do 

Soxt := 

I -.= 1 + 1 

for all 5 G S do 

for i = 1, . . . ,k^^^ do 

S' := S {copy partial solution S into 5'} 

S[ := ExtendPartialSolution(5',/,C) {see Algorithmic}- 

if solution S' is complete (that is, all tasks are assigned) then 

Bcompl '■— -Sconipl U {5 } 

else 

if ^ < m and S'l is different to the l-ih. work station of all other S € -Bcxt then 

-Boxt := -Bcxt U {S'] 
end if 
end if 
end for 
end for 

B ^SelectSolutions(i?cxt,fcbw) 
end while 

output: If -B^-on^pi 7^ the output is true, otherwise false 



of extensions /c^xt as input. Given a fixed cycle time C and m (the number of work stations) 
BS tries to find at least one feasible solution. As mentioned before, the crucial aspect of BS is 
the extension of partial solutions in several possible ways. At each step the algorithm extends 
each partial solution from i? in a maximum number of ways. More specifically, given a partial 
solution with I — 1 < m work stations already filled, an extension is generated by assigning 
a set of so-far unassigned tasks to the next work station Si such that the given cycle time 
C is not surpassed and the precedence constraints between the tasks are respected (see lines 
11-12 of Algorithm [T]) . The algorithm produces extensions in a (partially) probabilistic way 
rather than in the usual deterministic mannerjl Each generated extension (partial solution) 
is either stored in set .Bcompi in case it is a complete solution, or in set iJ^xt otherwise (see lines 
13-19 of Algorithm [1]) . However, a partial solution is only stored in set Bext if it uses at most 
m — 1 work stations, and if its l-th work station is different to the l-th work station of all 
partial solutions that are already in B. Finally, BS creates a new beam B by selecting up to 
/Cb„ solutions from set -B^xt of further extensible partial solutions (see line 22 of Algorithm [T]) . 
This is done in function SelectSolutions(-Bext;fcbw) on the basis of a lower bound function LB(-). 
In the following we describe in detail the extension of partial solutions and the working of 
function SelectSolutions(i?cxt)^bw)- 

^This is done for avoiding a choice-without-replacement process for which aU possible work station fillings 
must be generated beforehand. 
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Algorithm 2 Function ExtenclPartialSolution(5', /, C) of Algorithm [T] 
1: input: A partial solution S' , the index I of the work station to be filled, and the cycle 

time C 
2: S', := 



3 
4 
5 
6 
7 
8 
9 
10 
11 



r ■.= {ieT\i<^ Uj=i S'^,Pi'^ U=i S'j,ti + c_ < C} 

r ■= C 

while T' / do 

j :=ChooseTask(T', Crem) 

'-'rem • '-'rem 

S[ := SI U {j} 
end while 

output: Filled work station S'l 



Extending Partial Solutions. The generation of an extension of a partial solution S' with 
/ — 1 work stations already filled works as follows. Unassigned tasks are iteratively assigned 
to work station S'l until the sum of their processing times is such that no other task can be 
added to S'l without exceeding the given cycle time C. This procedure is pseudo-coded in 
Algorithm [2j At each step, T' denotes the set of so-far unassigned tasks that may be added 
to S'l without violating any constraints. The definition of this set of available tasks is given 
in line 3, respectively 8, of Algorithm [2j 

It remains to describe the implementation of function ChooseTask(T', c,^^) of Algorithm[2j 
For that purpose let us first define the following subset of T': 

T-' :={ieT'\ti + c,,^ = C} (8) 

This definition is such that T^""' contains all tasks that saturate, in terms of processing time, the 
Z-th work station Si. The choice of a task from T' is made on the basis of greedy information, 
that is, on the basis of values r/i > that are assigned to all tasks i € T' by a greedy function. 
The first action for choosing a task from T' consists in flipping a coin for deciding if the 
choice is made deterministically, or probabilistically. In case of a deterministic choice, there 
are two possibilities. First, if T*"^* 7^ 0, the best task from T'^''* is chosen, that is, the task with 
maximal greedy value among all tasks in T""**. Otherwise, we choose the task with maximal 
greedy value from T'. In case of a probabilistic decision, a task from T' is chosen on the basis 
of the following probability distribution: 

p(i):=^^ ,W€T' (9) 
jdT' 

For completing the description of function ChooseTask(T', c^m), we must describe the defini- 
tion of the greedy values ry,, Vi G T. In a first step a term 7j is defined as follows: 

ir.= ^i-{'-^)+n,-( '^^^'1' ,.1 ) ,V^G^ (10) 

\L J y maxi <j <„ I hnCj \ J 

Hereby, Sucf denotes the set of all tasks that can be reached from i in precedence graph G 
via a directed path. This definition combines two greedy heuristics that are often used in 
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the context of assembly line balancing problems. The first one concerns the task processing 
times and the second one concerns the size of Sucf . The influence of both heuristics can 
be adjusted via the setting of weights ki and K2- In order to be more flexible we decided to 
allow for both weights a value from [—1, 1]. This means that we consider for each heuristic 
potentially also its negation. Given the 7j-values, the greedy values rji are then derived as 
follows: 

^. .^T^^Trnin^ V^GT, (11) 

Tmax 

where %-nin, respectively j^^^, denote the minimum, respectively maximum, values of all 7^. 
Interestingly, for obtaining well- working greedy values, parameters ki and K2 have to be 
chosen in a problem-instance-dependent way. Tuning experiments are present in Section 14.21 



-'cxt • 



The Lower Bound Function. The new beam B is — at each step — chosen from 
This choice is implemented by function SelectSolutions(i?cxt)^bw) of Algorithm [2l First, the 
solutions in iJ^xt are ranked with respect to increasing lower bound values LB(-). Then, the 
minj/Cbw, |-B^xt|} highest ranked partial solutions from B^^^ are selected. Let us denote by 
T Q T the set of tasks that have not yet been assigned to work stations in partial solution 
S'. Then: 



LB(50 



C 



(12) 



Note that this lower bound is inspired by splitting-based bounds for the one-dimensional 
bin-backing problem. 

3.2 The Algorithmic Scheme 

The BS component outlined in the previous section is used by an iterative algorithmic 
scheme that is presented in Algorithm [3l Henceforth this algorithmic scheme is labelled 
iterated beam search (IBS). The first step consists in determining a starting cycle time C, 
which is computed in funcion DetermineStartingCycleTime() of Algorithm [3] as 



C := max < maxjgTjtj}, 



m 



(13) 



The algorithm works in two phases. In the first phase (see lines 3-11 of Algorithm [3D the 
algorithm tries to quickly find a first cycle time C for which a valid solution can be found. 
For this purpose BS is applied with the setting k^^ = 5 and k^^^ = 2. Note that this setting 
was chosen after tuning by hand. Moreover, note that the first phase only takes a fraction of a 
second of computation time. This holds for all instances considered in Section HI The second 
phase of the algorithm iteratively tries to find a valid solution for the next smaller cycle time. 
In this phase, the algorithm disposes over a certain time limit for each considered cycle time. 
Remember that the working of BS is partially probabilistic. Therefore, BS can repeatedly 
be applied to the same instance with potentially different outcomes. The first five percent of 
the above-mentioned time limit are spent by BS applications that use the setting fcbw := 10 
and /Ccxt := 5. This is done with the intention of not wasting too much computation time, if 
not necessary. However, if BS is not able to solve the given cycle time with this setting, the 
remaining 95% of the available time are spent by BS applications using the setting /c^w := 150 
and feext := 20. With this setting BS is considerably slower. However, the probability of 



7 



Algorithm 3 Iterative beam search (Ibs) for the SALBP-2 

1: input: an instance {T,G,m) 

2: C := DetermineStartingCycleTime() 

3: kf^^ := 5, fcgxt •= 2 

4: success := FALSE 

5: while not success do 

6: success := BeamSearch((r, G, m), C, /Cbw) ^ext) {original instance} 

7: if not success then 

8: success := BeamSearch{{T, ,m),C,k^^,ke^) {reverse instance} 

9: if not success then C := C + 1 end if 

10: end if 

11: end while 

12: C-C-1 

13: stop := FALSE 

14: while not stop do 

15: success := FALSE 

16: while time limit not reached and not success do 

17: if within 5% of time limit then fcb„ := 10, k^-^ := 5 else fcbw := 150, k^^ := 20 end 
if 

18: success := BeamSearch((T, G, m), C, febwj ^ext) {original instance} 

19: if not success then 

20: success := BeamSearch((r, C, m), C, fcbw, ^ext) {reverse instance} 

21: end if 

22: end while 

23: if success then C := C — 1 else stop := TRUE end if 

24: end while 

25: C :=C+1 

26: output: cycle time C 



finding feasible solutions is much higher than with the setting described before. The second 
phase of the algorithm ends when the time limit has passed without having found a feasible 
solution for the considered cycle time. 

4 Experimental Evaluation 

Ibs was implemented in ANSI C++, and GCC 3.4.0 was used for compiling the software. 
Experimental results were obtained on a PC with an AMD64X2 4400 processor and 4 Gb of 

memory. In the following we first describe the set of benchmark instances that wc used for 
the experimental evaluation. Subsequently we present the tuning process that we conducted 
and the results of the proposed algorithm. 

4.1 Benchmark Instances 

We used the usual set of 302 benchmark instances from the literature. They can be obtained — 

together with information about optimal and best-known solutions — from a website espe- 
cially dedicated to all kind of assembly line balancing problems maintained by Armin Scholl, 
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Table 1: Values of parameters ki 



and K2 for the final experiments. 



Graph 


Ki 




Graph 




K2 


Arcusl 


0.0 


1.0 


Lutz2 


-0.1 


0.9 


Arcus2 


-0.5 


0.9 


Lutz3 


0.0 


1.0 


Barthol2 


0.0 


1.0 


Mukherje 


0.0 


1.0 


Barthold 


0.0 


1.0 


Sawyer 


0.0 


1.0 


Buxey 


0.0 


1.0 


Scholl 


0.0 


1.0 


Gunther 


-0.4 


0.8 


Tonge 


-0.1 


0.2 


Hahn 


0.0 


1.0 


Warnecke 


0.0 


1.0 


Kilbridge 


0.0 


1.0 


Wee-Mag 


0.0 


1.0 


Lutzl 


-0.1 


0.9 









http://www.assembly-line-balanciiig.cie Each instance consists of a precedence graph 
G and a given number m of work stations. The benchmark set is composed of two subsets of 
instances, henceforth called Datasetl and Dataset2. Datasetl consists of 128 instances based 
on 9 different precedence graphs with a number of tasks between 29 to 111. Dataset2 is com- 
posed of 174 instances based on 8 different precedence graphs with a number of tasks varying 
from 53 to 297. 



4.2 Algorithm Tuning 

During preliminar experiments we realized that parameters febw and k^^t have a rather low 
impact on the final results of Ibs. In other words it is easy to find a reasonable setting 
for these parameters quite quickly. Their setting dynamically changes during a run of the 
algorithm as specified in Section 13.21 On the contrary, parameters ki and K2 (see Eq. [TOj) have 
a high impact on the algorithms' performance. Remember that ki is the weight of the greedy 
function concerning the task processing times, while K2 is the weight of the greedy function 
concerning the number of tasks that have to be processed after the task under consideration. 
As mentioned before, for both parameters we allowed values from [—1,1]. Instead of trying 
to find a good parameter setting for each instance, we decided for a tuning process based on 
precedence graphs, that is, we wanted to choose a single setting of ki and K2 for all instances 
concerning the same precedence graph. For that purpose we applied a specific version of Ibs 
for all combinations of ^1,^2 € { — 1-0, —0.9, . . . , 0.0, . . . , 0.9, 1.0} to all 302 instances. This 
makes a total of 441 different settings for each instance. The specific version of Ibs that 
we applied for the tuning process differs from Ibs as outlined in Algorithm [3] in that lines 
14-24 were replaced by a single, deterministic, application of beam search with k^,„ = 150 and 
fccxt = 20. This was done for the purpose of saving computation time. Based on the obtained 
results we chose the settings presented in Table [1] for the different precedence graphs. It 
is interesting to note that, apart from a few exceptions, the greedy heuristic based on task 
processing times does not seem necessary for obtaining good results. 

In Figure [T] tuning information is provided in graphical form for a few representative 
examples. The y-axis of the presented graphics varies over the different values of ki, while 
the X-axis ranges over the allowed values of K2. Note that each graphic consists of 441 squares 
representing the 441 different combinations of values for ki and K2- The gray level in which 
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each square is painted indicates the quahty of the algorithm when run with the corresponding 
parameter setting. In particular, black color denotes the worst setting, whereas white color 
indicates the best algorithm setting. In some cases such as (Buxey, m = 7) and (Tonge, 
m = 16), as shown in Figures 1(a) and 1(b), there is a wide range of good settings, which 
are basically all those with ni < 0. In other examples such as (Arcusl, m = 12) and (Scholl, 
m = 38) it is strictly required to set ki to and K2 to a positive value for obtaining good 



solutions; see Figures 1(c) and 1(d) Finally, the graphics shown in Figures 1(e) and 1(f) 



indicate that even for the same precedence graph a good parameter setting might depend 
strongly on the number of work stations. 



Table 2: Differences in algorithm performance when considering the best and the worst pa- 
rameter setting. 



Instance 


Best setting 


Worst setting 


Difference (%) 


(Buxey, m = 


7) 


47 


52 


10.64 


(Tonge, m = 


16) 


222 


265 


19.37 


(Arcusl, m = 


= 12) 


12599 


13767 


9.27 


(Scholl, m = 


38) 


1857 


2031 


9.37 


(Lutzl, m = 


9) 


1637 


1801 


10.02 


(Lutzl, m = 


10) 


1525 


1619 


6.16 



It is also interesting to quantify the differences in algorithm performance for different 
parameter settings. Table [2] shows for the six cases presented in Figure [1] the result of the 
algorithm with the best setting (column Best setting), the result of the algorithm with the 
worst setting (column Worst setting), and the difference (in percent) between these two 
settings. The results in Table [2] show that there are considerable differences in performance 
between the best and the worst algorithm setting. This underlines the importance of finding 
opportune values for ki and K2- 

4.3 Results 

Algorithm IBS was applied 20 times to all 302 instances. Herefore we used a computation 
time limit of 180 seconds for each cycle time, that is, Ibs was given maximally 180 seconds 
for finding a feasible solution for a given cycle time. In case of success, the algorithm has 
again 180 seconds for the next smaller cycle time, etc. Detailed results of IBS for all 302 
instances are given in Tables [7] and [8] that are to be found in Appendix A. The data is, in 
both tables, presented as follows. The first two columns provide the name of the precedence 
graph and the number of work stations (m). The third column provides the value of the 
optimal (respectively, best-known) solutions. In case a value is not proved to be optimal it 
is overlined. More in detail, in 15 out of 302 cases optimality has not been proved yet. The 
remaining five columns are reserved for the results of Ibs. The first of these five columns 
contains the value of the best solution found by Ibs over 20 runs. In case this value is 
presented with a gray background, a new best-known solution has been found. On the other 
side, if this value is marked by an asterisk, the obtained result does not reach the value of 
a best-known solution. In all other cases the values correspond to the values of best-known 
solutions. The second column provides the average over 20 runs, while the third column 
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-1.01 



-1.0 0.0 

(a) Buxcy. in = 




-1.0 0.0 

(c) Arcusl, m = 12 




-1.0 0.0 

(e) Lutzl, m = 9 




1.0 -1.0 0.0 

(b) Toug'c. Ill = 16 



0.0 



-1.0 















E 


■ 


1 


■ 



































1.0 -1.0 0.0 

(d) SchoU, m = 38 

1.0 



1.0 



0.0 



-1.0 



-1.0 0.0 

(f) Lutzl, m = 10 



1.0 



Figure 1: Tuning results presented in graphical form for six representative instances. The 
y-axis ranges over the values of ki, while the x-axis ranges over the values of K2- The gray 
levels of the squares indicate the quality of the algorithm when run with the corresponding 
setting: the lighter a square is painted, the better is the parameter setting. 
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contains the corresponding standard deviation. The fourth column gives the average time 
(in seconds) at which the best solution of a run was found, averaged over 20 runs. The fifth 
column provides the corresponding standard deviation. From the results presented in Tables [7| 
and [8] (see Appendix A) we can observe that Ibs obtains best-known (respectively, optimal) 
solutions in 276 out of 302 cases. Moreover, new best solutions are obtained in 6 cases. This 
is remarkable as — despite a considerable amount of ongoing research — in the last 14 years no 
improved solutions have been reported. Only in 20 cases (all concerning precedence graphs 
Arcusl, Arcus2, Scholl, and Warnecke) our algorithm was not able to find the best solutions 
known. However, in most of these cases the deviation from the best-known solution is no 
more than one unit of cycle time. 



Table 3: Results of Ibs in comparison to the best (TabuSearch), respectively most recent 
(DE_RKS and PNA-for), methods from the literature. 





DE_RKS 


PNA-for 


TabuSearch 


Iterative Beam Search (Ibs) 


#opt 


n/g 


39 


168 


282 


mrd (%) 


2.64 


2.73 


0.40 


0.0029 


time (s) 


10.74 


407.28 


84.90 


31.61 


Note: n/g meant 


not given 









In addition to Tables [7] and [8] the results of IBS are presented in a summarized way in 
Table [3l in comparison to three other algorithms. TabuSearch [22j, even though already 
published in 1996, still counts as the current state-of-the-art heuristic method for SALBP-2. 
DE_RKS is the best version of a differential evolution (DE) algorithm proposed in [I4j, and 
PNA-FOR is a Petri net-based heuristic published in |llj . The last two methods are, to our 
knowledge, the most recently published heuristic methods for SALBP-2. Three measure are 
used in Table [3] for the comparison of Ibs with these three algorithms. The row labelled #opt 
provides the number of best-known solutions found by each method (over 302). Moreover, the 
row labelled mrd (%) gives the mean relative deviation (in percent) of the results obtained 
by the four algorithms from the best-known solutions for all 302 instances. Finally, row time 
contains the average computation time of the algorithms for all 302 instances. Concerning 
the quality of the results, we can conclude that Ibs clearly outperforms its competitors. For 
the correct interpretation of the computation times it has to be taken into account that the 
four algorithms were executed on computers with very different processor speeds. While 
TabuSearch was executed on a 80486 DX2-66 processor, PNA-FOR was run on an Athlon 
XP 2000+ processor with 1.67 GHz, and DE_rks was run on a Pentium IV processor with 
1.7 GHz. This means that TabuSearch was run by far on the slowest machine, Ibs by 
far on the fastest machine, and PNA-FOR and DE_rks on comparable machines. Given 
the computation times in Table [3] we can savely conclude that TabuSearch is the fastest 
algorithm, and PNA-FOR is the slowest one. However, note that assembly line balancing is, 
in most cases, not a time-critical application. In other words, for most practical purposes it 
does not matter if an algorithm takes 1 minute or 6 hours of computation time. 

In the following we present the results of Ibs in comparison to DE_rks and PNA-FOR in 
the same way as done in [H] and [11]. In these works, results were presented as averages over 
instances based on the same precedence graph, and also averaged over Datasetl and Dataset2. 
The quality of the results is given in terms of the mean relative deviation (in percent) from 
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Table 4: Results of Ibs in comparison to DE_rks and PNA-FOR for the 128 instances of 
Datasetl (averaged over precedence graphs). 



Graph 


mrd 


DE_RKS 

(%) time (s) 


PNA-FOR 
mrd (%) 


Iterative Beam Search (Ibs) 
mrd (%) time (s) 


Buxey 


1.16 


0.80 


3.07 


0.0 


0.06 


Sawyer 


2.27 


1.64 


4.00 


0.0 


0.14 


Lutzl 


0.32 


0.88 


4.09 


0.0 


0.55 


Gunther 


0.14 


1.08 


1.27 


0.0 


0.02 


Kilbridge 


0.66 


1.43 


2.19 


0.0 


0.0067 


Tonge 


1.88 


3.71 


2.53 


0.0 


7.73 


Arcusl 


0.99 


5.29 


2.47 


0.0287 


152.34 


Lutz2 


3.08 


1.00 


2.99 


0.0 


0.035 


Arcus2 


4.96 


19.02 


2.06 


0.0066 


135.83 


Average: 


1.72 


3.87 


2.57 


0.0058 


51.76 



the best-known solutions. Tables H] and [5] clearly show that Ibs is largely superior to both 
competitor algorithms. 

4.4 Results of a High-Performance Version 

In an attempt to further improve on the results of our algorithm we decided to apply a 
high-performance version of IBS to all problem instances for which the optimal solution is 
unknown and, additionally, to all instances where IBS — with the settings as outlined in the 
previous section — was not able to find a best-known solution. This high-performance version 
is obtained as follows. First, 1800 seconds are used as a time limit for each cycle time. 
Second, in line 17 of Algorithm [3] only 1% of the time limit is used (instead of 5%). Third, 
for each application of beam search in lines 18 and 20 of Algorithm [3] the beam width /c^w 
is randomly chosen from [150, 250] and the number of extensions is randomly chosen from 
[20, 40]. Moreover, with a probability of 0.5 the heuristic information is — for each application 
of beam search — calculated using the weight values as outlined in Table [1] Otherwise, the 
weight values are chosen randomly from [—1,1]. With these modifications we applied IBS 
exactly once to all the instances of Table [H The results of the algorithm are given in column 
Result. Indeed, the number of instances for which a best-known solution can not be found 
before is reduced from 20 to 10 instances. Moreover, the algorithm is now able to find 
new best-known solutions in 9 (instead of only 6) cases. Summarizing, this amounts to 283 
best-known solutions found and 9 new best-known solutions obtained. In one case, (Scholl, 
m = 49), the new best-known solution is provenly optimal, as its value coincides with the 
best known lower bound. 

5 Conclusions and Future Work 

In this work we have proposed an iterative beam search algorithm for the simple assembly 
line balancing problem with a fixed number of work stations, SALBP-2. The experimental 
evalution of the algorithm has shown that it is currently a state-of-the-art method for this 
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Table 5: Results of Ibs in comparison to DE_rks and PNA-FOR for the 174 instances of 
Dataset2 (averaged over precedence graphs). 



Graph 


mrd 


DE_RKS 

(%) time (s) 


PNA-FOR 
mrd (%) 


Iterative Beam Search (Ibs) 
mrd (%) time (s) 


Hahn 


0.0 


1.00 


2.52 


0.0 


0.065 


Warnecke 


3.74 


3.53 


5.57 


0.0579 


1.54 


Wee-Mag 


1.23 


3.68 


1.56 


0.0 


0.65 


Lutz3 


1.68 


5.62 


2.59 


0.0 


0.61 


Mukherje 


n/a 


n/a 


1.04 


0.0 


2.52 


Barthold 


0.26 


19.47 


1.02 


0.0 


0.079 


Barthol2 


6.85 


33.08 


3.97 


0.0 


0.96 


Scholl 


9.51 


43.95 


3.21 


0.0028 


98.69 


Average: 


3.32 


15.79 


2.85 


0.00071 


16.79 



problem. Appart from producing optimal, respectively best-known, solutions in 283 out of 
302 test cases, our algorithm generated new best-known solutions in further 9 test cases. 
Encouraged by the results for the SALBP-1 version of the problem (as published in |3j) and 
the results obtained in this paper for the SALBP-2 we intent to apply similar algorithms 
based on beam search to other assembly line balancing problems. 
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Table 7: Detailed results of Ibs for 302 test instances (Part A). 
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best 


avg 


std 


time (s) 


std 








best 


avg 


std 


time (s) 


std 


Arcusl (83) 


3 


25236 


25236 


25236.00 


(0.00) 


0.56 


(0.29) 


Barthold (148) 


9 


626 


626 


626.00 


(0.00) 


0.13 


(0.03) 




4 


18927 


18927 


18927.00 


(0.00) 


6.16 


(3.92) 




10 


564 


564 


564.00 


(0.00) 


0.10 


(0.02) 




B 


15142 


15142 


15142.00 


(0.00) 


80.08 


(16.69) 




11 


513 


513 


513.00 


(0.00) 


0.09 


(0.00) 




6 


12620 


12620 


12620.00 


(0.00) 


15.52 


(5.76) 




12 


470 


470 


470.00 


(0.00) 


0.09 


(0.00) 




7 


10826 


10826 


10826.00 


(0.00) 


59.80 


(19.89) 




13 
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434 


434.00 


(0.00) 


0.16 


(0.02) 




8 


9554 


•9555 


9555.55 


(0.51) 


120.54 


(63.16) 




14 


403 


403 


403.00 


(0.00) 


0.14 


(0.01) 




9 


8499 


8499 


8500.70 


(0.80) 


130.73 


(117.23) 




15 


383 


383 


383.00 


(0.00) 


0.05 


(0.00) 




10 


7580 


7580 


7580.95 


(0.60) 


311.89 


(109.87) 


Buxey (29) 


7 


47 


47 


47.00 


(0.00) 


0.46 


(0.51) 




11 


7084 


•7086 


7086.70 


(0.47) 


71.13 
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41 


41 


41.00 
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6414.10 
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37 


37 
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0.01 


(0.00) 




13 
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5864.00 


(0.00) 


129.18 


(12.99) 
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34 


34 


34.00 
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0.00 


(0.00) 




14 


5441 


5441 


5441.00 


(0.00) 


2.02 


(0.41) 
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0.00 
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(0.00) 




16 


4850 


4850 


4850.00 
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(0.00) 


0.47 


(0.05) 




14 


35 


35 


35.00 


(0.00) 


0.03 


(0.00) 




35 


121 


121 


121.00 


(0.00) 


0.94 


(0.15) 




15 


33 


33 


33.00 


(0.00) 


0.03 


(0.00) 




36 


118 


118 


118.00 


(0.00) 


0.87 


(0.11) 




16 


31 


31 


31.00 


(0.00) 


0.02 


(0.00) 




37 


115 


115 


115.00 


(0.00) 


0.92 


(0.19) 




17 


29 


29 


29.00 


(0.00) 


0.04 


(0.00) 




38 


112 


112 


112.00 


(0.00) 


0.76 


(0.06) 




18 


28 


28 


28.00 


(0.00) 


0.04 


(0.00) 




39 


109 


109 


109.00 


(0.00) 


0.70 


(0.13) 




19 


26 


26 


26.00 


(0.00) 


0.04 


(0.01) 




40 


106 


106 


106.00 


(0.00) 


0.91 


(0.32) 




20 


25 


25 


25.00 


(0.00) 


0.04 


(0.00) 




41 


104 


104 


104.00 


(0.00) 


0.72 


(0.05) 




21 


24 


24 


24.00 


(0.00) 


0.02 


(0.00) 




42 


101 


101 


101.00 


(0.00) 


2.77 


(5.43) 




22 


23 


23 


23.00 


(0.00) 


0.04 


(0.00) 




43 


99 


99 


99.00 


(0.00) 


0.84 


(0.23) 




23 


22 


22 


22.00 


(0.00) 


0.05 


(0.00) 




44 


97 


97 


97.00 


(0.00) 


0.59 


(0.15) 




24 


21 


21 


21.00 


(0.00) 


0.05 


(0.00) 




45 


95 


95 


95.00 


(0.00) 


0.76 


(0.25) 




25 


20 


20 


20.00 


(0.00) 


0.02 


(0.00) 




46 


93 


93 


93.00 


(0.00) 


0.68 


(0.29) 




26 


19 


19 


19.00 


(0.00) 


0.05 


(0.02) 




47 


91 


91 


91.00 


(0.00) 


0.67 


{(),1.">) 




27 


19 


19 


19.00 


(0.00) 


0.03 


(0.00) 




48 


89 


89 


89.00 


(0.00) 


1.17 


(0.47) 




28 


18 


18 


18.00 


(0.00) 


0.05 


(0.00) 




49 


87 


87 


87.00 


(0.00) 


0.83 


(0.26) 


Lutz3 (89) 


3 


548 


548 


548.00 


(0.00) 


0.01 


(0.00) 




50 


85 


85 


85.90 


(0.31) 


3.28 


(10.33) 




4 


411 


411 


411.00 


(0.00) 


0.03 


(0.00) 




51 


84 


84 


84.00 


(0.00) 


2.60 


(2.09) 




5 


329 


329 


329.00 


(0.00) 


0.01 


(0.00) 


Barthold (148) 


3 


1878 


1878 


1878.00 


(0.00) 


0.01 


(0.00) 




6 


275 


275 


275.00 


(0.00) 


0.01 


(0.00) 




4 


1409 


1409 


1409.00 


(0.00) 


0.02 


(0.00) 




7 


236 


236 


236.00 


(0.00) 


0.02 


(0.00) 




5 


1127 


1127 


1127.00 


(0.00) 


0.02 


(0.00) 




8 


207 


207 


207.00 


(0.00) 


0.04 


(0.00) 




6 


939 


939 


939.00 


(0.00) 


0.06 


(0.01) 




9 


184 


184 


184.00 


(0.00) 


5.54 


(4.14) 




7 


805 


805 


805.00 


(0.00) 


0.11 


(0.00) 




10 


165 


165 


165.00 


(0.00) 


0.09 


(0.01) 




8 


705 


705 


705.00 


(0.00) 


0.04 


(0.01) 




11 


151 


151 


151.00 


(0.00) 


0.06 


(0.00) 
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Table 8: Detailed results of Ibs for 302 test instances (Part B). 



Graph 


m 


bks 




Iterative Beam Search (IBS 




Graph 




m. 


bks 




Iterative Beam Search (IBS) 








best 


avg 


std 


time (s) 


std 










best 


avg 


std 


time (s) 


std 


Lutz3 (89) 


12 


138 


138 


138.00 


(0.00) 


0.02 


(0.00) 


Tonge 


(70) 


6 


,585 


,585 


.585.00 


(0.00) 


0,03 


(0,02) 




13 


128 


128 


128.00 


(0.00) 


0.05 


(0.03) 






7 


502 


502 


,502.00 


(0.00) 


0,06 


(0,01) 




14 


118 


118 


118.00 


(0.00) 


0.67 


(0.52) 






8 


439 


439 


439.00 


(0.00) 


0,04 


(0,01) 




15 


110 


110 


110.00 


(0.00) 


0.80 


(0.66) 






9 


391 


391 


391.00 


(0,00) 


0,05 


(0,00) 




16 


105 


105 


105.00 


(0.00) 


0.34 


(0.30) 






10 


352 


352 


352.00 


(0,00) 


0,06 


(0,01) 




17 


98 


98 


98.00 


(0.00) 


0.19 


(0.05) 






11 


320 


320 


320.00 


(0,00) 


0,06 


(0,02) 




18 


93 


93 


93.00 


(0.00) 


0.16 


(0.05) 






12 


294 


294 


294.00 


(0,00) 


0,09 


(0,01) 




19 


89 


89 


89.00 


(0.00) 


0.05 


(0.03) 






13 


271 


271 


271.00 


(0,00) 


0,06 


(0,01) 




20 


85 


85 


85.00 


(0.00) 


0.14 


(0.05) 






14 


251 


251 


251.80 


(0,41) 


20,25 


(45,26) 




21 


80 


80 


80.00 


(0.00) 


0.11 


(0.02) 






15 


235 


235 


235.00 


(0,00) 


0,35 


(0,16) 




22 


76 


76 


76.00 


((>.()()) 


4.20 


(3.30) 






16 


221 


221 


221.00 


(0,00) 


0,97 


(0,93) 




23 


74 


74 


74.00 


(0.00) 


0.29 


(0.18) 






17 


208 


208 


208.00 


(0,00) 


1,02 


(0,95) 


Mukherje (94) 


3 


1403 


1403 


1403.00 


(0.00) 


0.01 


(0.00) 






18 


196 


196 


196.15 


(0.37) 


46,76 


(39,05) 




4 


1052 


1052 


1052.00 


(0.00) 


0.10 


(0.00) 






19 


186 


186 


186.35 


(0.49) 


87,09 


(64,62) 




5 


844 


844 


844.00 


(0.00) 


0.02 


(0.00) 






20 


177 


177 


177.90 


(0.31) 


9,46 


(39,39) 




6 


704 


704 


704.00 


(0.00) 


0.02 


(0.00) 






21 


170 


170 


170.00 


(0.00) 


8,71 


(9,43) 




7 


621 


621 


621.00 


(0.00) 


0.12 


(0.00) 






22 


162 


162 


162.00 


(0.00) 


2,52 


(2,45) 




8 


532 


532 


532.00 


(0.00) 


0.09 


(0.01) 






23 


156 


156 


156.00 


(0.00) 


0,11 


(0,04) 




9 


471 


471 


471.00 


(0.00) 


0.09 


(0.01) 






24 


156 


156 


156.00 


(0.00) 


0,02 


(0,00) 




10 


424 


424 


424.00 


(0.00) 


0.02 


(0.00) 






25 


156 


156 


156.00 


(0,00) 


0,02 


(0,00) 




11 


391 


391 


391.00 


(0.00) 


0.10 


(0.01) 


Warnecke (58) 


3 


516 


516 


516.00 


(0,00) 


0,02 


(0,00) 




12 


358 


358 


358.00 


(0.00) 


0.05 


(0.01) 






4 


387 


387 


,387,00 


(0,00) 


0,08 


(0,08) 




13 


325 


325 


325.00 


(0.00) 


0.17 


(0.06) 






5 


310 


310 


310,00 


(0.00) 


0,03 


(0,00) 




14 


311 


311 


311.00 


(0.00) 


0.08 


(0.01) 






6 


258 


258 


2,58,00 


(0.00) 


0,10 


(0,06) 




15 


288 


288 


288.00 


(0.00) 


0.09 


(0.02) 






7 


222 


222 


222,00 


(0.00) 


0,03 


(0,00) 




16 


268 


268 


268.00 


(0.00) 


0.20 


(0.02) 






8 


194 


194 


194,00 


(0.00) 


0,05 


(0,01) 




17 


251 


251 


251.00 


(0.00) 


0.45 


(0.05) 






9 


172 


172 


172,00 


(0.00) 


2,40 


(2,28) 




18 


239 


239 


239.00 


(0.00) 


0.24 


(0.02) 






10 


155 


155 


155,00 


(0.00) 


0,05 


(0,01) 




19 


226 


226 


226.00 


(0.00) 


0.03 


(0.01) 






11 


142 


142 


142,00 


(0.00) 


0,04 


(0,00) 




20 


220 


220 


220.05 


(0.22) 


57.07 


(43.23) 






12 


130 


130 


130,00 


(0.00) 


0,07 


(0,01) 




21 


208 


208 


208.00 


(0.00) 


0.12 


(0.00) 






13 


120 


120 


120,00 


(0.00) 


0,08 


(0,01) 




22 


200 


200 


200.00 


(0.00) 


0.13 


(0.06) 






14 


111 


111 


111,00 


(0.00) 


0,29 


(0,24) 




23 


189 


189 


189.00 


(O.OO) 


0.12 


(0.02) 






15 


104 


104 


104,00 


(0.00) 


0,27 


(0,13) 




24 


179 


179 


179.00 


(().()()) 


0.48 


(0,41) 






16 


98 


98 


98,00 


(0.00) 


0,08 


(0,01) 




25 


172 


172 


172.00 


(0.00) 


0.43 


(0.18) 






17 


92 


92 


92,00 


(0.00) 


0,26 


(0,14) 




26 


171 


171 


171.00 


(0.00) 


0.27 


(0,03) 






18 


87 


87 


87,00 


(0,00) 


0,27 


(0,14) 


Sawyer (30) 


7 


47 


47 


47.00 


(0.00) 


0.14 


(0,11) 






19 


84 


84 


84,00 


(0,00) 


0,10 


(0,06) 




8 


41 


41 


41.00 


(0.00) 


0.06 


(0,03) 






20 


79 


79 


79,00 


(0,00) 


0,13 


(0,03) 




9 


37 


37 


37.00 


(0.00) 


0.03 


(0,01) 






21 


76 


76 


76,00 


(0,00) 


0,08 


(0,04) 




10 


34 


34 


34.00 


(0.00) 


0.02 


(0,00) 






22 


73 


73 


73,00 


(0,00) 


0,53 


(0,36) 




11 


31 


31 


31.00 


(0.00) 


0.82 


(0,85) 






23 


69 


69 


69,00 


(0,00) 


18,10 


(21,25) 




12 


28 


28 


28.00 


(0.00) 


0.03 


(0,01) 






24 


66 


66 


66,90 


(0,31) 


2,88 


(9,10) 




13 


26 


26 


26.00 


(0.00) 


0.02 


(0,00) 






25 


64 


•65 


65,00 


(0,00) 


12,30 


(11,87) 




14 


25 


25 


25.00 


(0.00) 


0.02 


(0,00) 






26 


64 


64 


64,00 


(0,00) 


0,05 


(0,01) 


SchoU (297) 


25 


2787 


2787 


2787.00 


(0.00) 


37.11 


(6,11) 






27 


60 


60 


60,00 


(0,00) 


1,61 


(0,73) 




26 


2680 


2680 


2680.00 


(0.00) 


32.56 


(4,71) 






28 


59 


59 


,59,00 


(0,00) 


0,16 


(0,06) 




27 


2580 


2580 


2580.00 


(0.00) 


58.06 


(8,65) 






29 


56 


56 


56,00 


(0,00) 


1,43 


(1,15) 




28 


2488 


2488 


2488.00 


(0.00) 


53.78 


(12,43) 


Wee-M 


ag (75) 


3 


500 


500 


500,00 


(0.00) 


0,01 


(0,00) 




29 


2402 


2402 


2402.00 


(0.00) 


45.57 


(9,93) 






4 


375 


375 


375,00 


(0.00) 


0,03 


(0,00) 




30 


2322 


2322 


2322.30 


(0.47) 


96.86 


(53,36) 






5 


300 


300 


300,00 


(0.00) 


0,02 


(0,00) 




31 


2247 


•2248 


2248.00 


(0.00) 


51.23 


(4,32) 






6 


250 


250 


250,00 


(0.00) 


0,03 


(0,00) 




32 


2177 


2177 


2177.50 


(0.51) 


91.50 


(63,17) 






7 


215 


215 


215,00 


(0,00) 


0,01 


(0,00) 




33 


2111 


2111 


2111.80 


(0.41) 


57.19 


(39,95) 






8 


188 


188 


188,00 


(0,00) 


0,04 


(0,01) 




34 


2049 


2049 


2049.20 


(0.41) 


117.47 


(53,03) 






9 


167 


167 


167,00 


(0,00) 


0,01 


(0,00) 




35 


1991 


1991 


1991.00 


(0.00) 


95.54 


(17,19) 






10 


150 


150 


1,50,00 


(0,00) 


0,08 


(0,02) 




36 


1935 


•1936 


1936.00 


(0.00) 


61.22 


(8,78) 






11 


137 


137 


137,00 


(0,00) 


0,06 


(0,01) 




37 


1883 


1883 


1883.05 


(0.22) 


140.57 


(34,46) 






12 


125 


125 


125,00 


(0,00) 


0,10 


(0,02) 




38 


1834 


1834 


1834.00 


(0.00) 


125.92 


(34,53) 






13 


116 


116 


116,00 


(0,00) 


0,02 


(0,00) 




39 


1787 


1787 


1787.35 


(0.49) 


192.12 


(55,85) 






14 


108 


108 


108,00 


(0,00) 


0,54 


(0,60) 




40 


1742 


1742 


1742.00 


(0.00) 


153.02 


(29,66) 






15 


100 


100 


100,00 


(0,00) 


0,24 


(0,19) 




41 


1700 


1700 


1700.00 


(0.00) 


81.51 


(18,06) 






16 


94 


94 


94,00 


(0,00) 


0,05 


(0,01) 




42 


1659 


•1660 


1660.00 


(0.00) 


146.58 


(32,11) 






17 


89 


89 


89,00 


(0,00) 


0,42 


(0,28) 




43 


1621 


1621 


1621.00 


(0.00) 


149.08 


(21,75) 






18 


87 


87 


87,00 


(0,00) 


0,07 


(0,01) 




44 


1584 


1584 


1584.15 


(0.37) 


127.52 


(54,06) 






19 


85 


85 


85,00 


(0,00) 


0,08 


(0,02) 




45 


1549 


1549 


1549.00 


(O.OO) 


1.50.65 


(32,59) 






20 


77 


77 


77,00 


(0,00) 


0,13 


(0,04) 




46 


1515 


•1516 


1516.00 


(0.00) 


121.28 


(20,06) 






21 


72 


72 


72,00 


(0,00) 


0,08 


(0,01) 




47 


1483 


•1484 


1484.00 


(0.00) 


165.44 


(28,64) 






22 


69 


69 


69,00 


(0,00) 


0,06 


(0,01) 




48 


1452 


1452 


1452.70 


(0.47) 


151.98 


(60,76) 






23 


67 


67 


67,00 


(0,00) 


4,34 


(3,58) 




49 


1427 




1424.00 


(0.00) 


113.46 


(22,47) 






24 


66 


66 


66,00 


(0,00) 


0,25 


(0,14) 




50 


1394 


1394 


1,394.85 


(0.37) 


143.63 


(.53,48) 






25 


65 


65 


65,00 


(0,00) 


11,36 


(26,58) 




51 


1386 


1386 


1386.00 


(O.OO) 


2.23 


(0,09) 






26 


65 


65 


65,00 


(0,00) 


0,05 


(0,01) 




52 


1386 


1386 


1386.00 


(0.00) 


0.17 


(0,01) 






27 


65 


65 


65.00 


(0,00) 


0,02 


(0,00) 


Tonge (70) 


3 


1170 


1170 


1170.00 


(0.00) 


0.01 


(0,00) 






28 


64 


64 


64.00 


(0,00) 


0,02 


(0,00) 




4 


878 


878 


878.00 


(0.00) 


0.02 


(0,00) 






29 


63 


63 


63.00 


(0,00) 


0,02 


(0,00) 




5 


702 


702 


702.00 


(0.00) 


0.05 


(0,02) 






30 


56 


56 


56.00 


(0,00) 


0,07 


(0,00) 
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